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We investigate density fluctuations in a coherent ensemble of interacting fermionic atoms. Adapt- 
ing the concept of full counting statistics, well-known from quantum optics and mesoscopic electron 
transport, we study second-order as well as higher-order correlators of density fluctuations. Using 
the mean-field BCS state to describe the whole interval between the BCS limit and the BEC limit, 
we obtain an exact expression for the cumulant-generating function of the density fluctuations of 
an atomic cloud. In the two-dimensional case, we obtain a closed analytical expression. Poissonian 
fluctuations of a molecular condensate on the BEC side are strongly suppressed on the BCS side. 
The size of the fluctuations in the BCS limit is a direct measure of the pairing potential. We also 
discuss the BEC-BCS crossover of the third cumulant and the temperature dependence of the second 
cumulant. 



I. INTRODUCTION 

Following the successful creation of Bose- Einstein con- 
densates (BECs) in ultracold atomic clouds [l|, re- 
cently ultracold Fermionic clouds have been produced 
@; H, BL IS ■ This has attracted a lot of attention both 
theoretically and experimentally, especially due to the 
ability to tune the mutual interaction between atoms via 
a Fano-Feshbach resonance. The unique opportunity to 
study the crossover from weak attractive to strong at- 
tractive interactions in one and the same system makes 
this interesting from a fundamental many-body point of 
view. 

Theoretically, fermionic systems with weak attractive 
interaction are superfluids and as such described by the 
Bardeen-Cooper-Schrieffer (BCS) theory Q. This theory 
can also describe the limit of stronger attractive interac- 
tion [13], in which a BEC of molecules is formed. In 
the crossover regime, the long-range nature of the inter- 
action makes the BCS theory less accurate [Tl| . 

Recently, several experiments studied the strongly in- 
teracting BEC-BCS crossover regime using spin mixtures 
of ultracold fermionic gases (see 12[ for a recent re- 
view). Measurements of the interaction energy of an 
ultracold fermionic gas near a Feshbach resonance were 
made, studying the impact of the interaction on the time 
of flight expansion [131] . Experimental investigation of 
collective excitations showed a strong dependence on the 
cou plin g strength [3, EH . Moreover, the condensation 
[iH, and the spatial correlations [l8| of the fermionic 
atom pairs were observed in the full crossover regime. 
The pairing gap was measured directly via a spectro- 
scopic technique in the whole crossover region UM- Re- 
markably, the gap values are in relatively good agreement 
with the BCS model in the whole region. 

The use of noise correlations to probe the many-body 
states of ultracold atoms was proposed in [2(| (see also 
pl|). Correlation measurements can be applied to de- 
tect phase coherence in mesoscopic superpositions [22j . 



The density and spin structure factor for the BEC-BCS 
transition was calculated (23[. Interferometric measure- 
ment schemes of the spatial pairing order have been pro- 
posed based on the atom counting statistics in the output 
channels [HI • Pairin g flu ctuations of trapped Fermi gases 
have been studied in 25| . The number statistics of Fermi 
and Bose gases has also been investigated in (2|| H3| • In 
recent experiments [28] the spatial structure of an atomic 
cloud has been directly observed (without the expansion 
used in most other experiments). This makes it possible 
to determine the density fluctuations either by repeat- 
ing the experiment many times or by taking densities at 
different positions in a homogeneous system to extract 
the statistics. Atomic shot noise has been experimen- 
tally investigated both in bosonic and fermionic systems 

[II Hil [Hi 11, [Mill- 



In this article we propose to use the full counting statis- 
tics of density fluctuations as a tool to gain access to 
the many-body nature of the ground state of a fermionic 
cloud in the BEC-BCS crossover regime (for other work 
on full counting statistics in ultracold atomic gases see 
[36[ on the experimental and [13, [H, [H, H(| on the the- 
oretical side). Our main result is a general expression 
for the particle number statistics of the mean-field BCS 
wave function. In the limiting cases, the statistics allows 
a straightforward interpretation. Deep in the molecular 
BEC limit the statistics is Poissonian, i.e. that of in- 
dependent pairs of atoms. In the opposite limit, on the 
BCS side of the crossover, the fluctuations are strongly 
suppressed and reflect the particle-hole symmetry. The 
statistics in the crossover regime differs strongly from 
both the BCS and BEC limits and will be discussed based 
on several numerical results. 
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II. COUNTING STATISTICS OF DENSITY 
FLUCTUATIONS 

We consider an atomic cloud with spatial distribution 
n(x) and divide the system into bins, see Fig. [TJ We 
are interested in the probability P{N) — (5(N — N)) 
to find N atoms in a "bin" of the system, see Fig. [T] 
Here, N = J v h(x) is the bin number operator, Vbin is 
the bin volume, and n(x) denotes the density operator. 
The bin volume is assumed to be much smaller than the 
volume of the full system, but still much larger than the 
mean interparticle distance cubed. Hence, a single bin 
can be considered as a grand-canonical system with the 
surrounding atomic cloud serving as the reservoir. 

In practice, it is more convenient to study the Fourier 
transform of P(N), which is related to the cumulant gen- 
erating function (CGF) S(x) via 

e ~ S(x) = 2 e txN P(N) = {e lxfl ) . (1) 

JV 

The cumulants C n are defined in a standard way as 
S(x) = — J2 n C n (ix) n / n * an d can be used to charac- 
terize the full counting statistics (FCS). We recall that 
C\ = N is the average number of atoms in a bin, and 
C*2 = ((N — N) 2 } measures the width of the number 
distribution. The third cumulant is proportional to the 
skewness and, therefore, a measure of the asymmetry of 
the distribution function. As we will see later, the third 
cumulant in a fermionic system is related to deviations 
from particle-hole symmetry. 




FIG. 1: (Color online) (a) Sketch of a typical atomic number 
density. The measured observable is the atom number TV 
given by n(x) integrated over the bin volume (the bins are 
indicated by the dashed lines), (b) Histogram of probabilities 
to find N particles in a bin. 



III. NONINTERACTING FERMIONS AND 
BOSONS 

We start by recalling some properties of noninteracting 
Fermi and Bose gases at a given temperature T. Here, 
individual atoms are independent and we can obtain the 
statistics as 

S(x)=T£ln[l±/±(*)(e ix -l)] - ( 2 ) 

k 

where f±(k) — [exp((efc — fij/ksT) ± l]" 1 is the Fermi 
(Bose) occupation function. The chemical potential is 



determined by the average atom number, i. e., C\ — 
N = /±- Obviously, the number statistics of fcrmions 
and bosons differs drastically in the degenerate regime. 
In the nondegenerate regime f± is small and we find 
the Poissonian statistics of classical particles S(x) — 
—N(exp(ix) — 1), both for fermions and bosons. 
For degenerate fermions ksT <C £f the statistics is 

S( x ) = -i X N - (Dk B T/Ae F )N X 2 ■ (3) 

Here and in the following, S(x) is defined in the interval 
[— 7T, 7r] and extended periodically. Thus, particle-number 
fluctuations are suppressed by T j tp in comparison to the 
classical case [4l[ . Remarkably, the statistics is Gaussian 
and consequently all cumulants C n for n > 3 vanish. This 
behavior resembles fermions in a ID wire [42| and can be 
interpreted as a consequence of anti-bunching. Note, that 
the Gaussian nature of the statistics is a consequence 
of particle-hole symmetry and is therefore strictly lim- 
ited to the degenerate regime. Higher-order corrections 
in ksT/eF will introduce deviations from the Gaussian 
limit, and lead e.g. to the appearance of higher-order 
odd cumulants, which are directly related to deviations 
from the perfect particle-hole symmetry. 

In contrast to that, we obtain quite a different behav- 
ior for free bosons. Approaching the degeneracy tem- 
perature T<? EC = 27r?i 2 n 2 / 3 /rnC(3/2) 2/3 [ll from above, 
the fluctuations are enhanced due to the large factor 
f-(k)(f-(k) + 1). In the condensed regime, the oc- 
cupation of the ground state becomes macroscopically 
larg e and the grand-canonical approach is no longer valid 
HIT EH. This can be seen if we take the limit T = of 
Eq. ©, leading to S(x) = In [l - N(e lx - 1)], which 
corresponds to a negative binomial distribution and the 
fluctuations therefore diverge according to C n ~ N n . As 
in the grand-canonical ensemble the chemical potential is 
[X = below the critical temperature T < Tjl? EC , an ar- 
bitrarily large number of bosons can be transferred from 
the reservoir into the bin, leading to unphysically large 
fluctuations. In order to find the correct fluctuations in 
this case, we calculate the FCS from Eq. ([T]) explicitly. 
We divide the boson operator into au = + , where 

b k = [ d 3 r e lfcr *(r) , c k = [ d 3 r e ikr ^{r) , (4) 

where V is the volume V without the volume Vbin of 
the bin. We consider a fully condensed state of A^ot 
non-interacting bosons \ip) — (aj)^* * |vac) = (fej + 
Co) Ntot | vac). Using Eq. (fTJ) for the bin number opera- 
tor N = b k bk we obtain a binomial statistics for a 
non-interacting bosonic gas in a bin at T = 

S(x) = -JVtotln[l + ^(e^-l)]. (5) 

For bin volumes Vb- m /V <C 1 the particle- number statis- 
tics becomes Poissonian, i.e. S(x) ~ —N(e lx — 1), where 
= N tot V hin /V. 



IV. BCS GROUND STATE 

We now turn to a Fermi gas with an attractive interac- 
tion parameterized by a scattering length a. The Hamil- 
tonian for the full system is number conserving and as 
such shows no number fluctuations at all. Here we con- 
sider a bin, i.e. a small subsystem, which we assume to 
be described by the non-number conserving BCS Hamil- 
tonian A straightforward ansatz for a non-number 
conserving many-body state, which takes the pairing in- 
teraction into account, is the BCS wave function, which 
will be used in the following and is known to correctly 
describe both the BCS and the BEC limit We 
will later show that the approach also reproduces the cor- 
rect counting statistics in the two limits, and therefore 
prefer to use this transparent, almost analytical approach 
instead of more complex canonical approaches. The BCS 
wave function is given by 

|BCS>=n(«fc + ^4A i )|0> . (6) 

k 

The variational procedure yields v 2 = 1 — u| = 
(1 - (e k - n)/E k ) 1% where E\ = (e k - ^) 2 + A 2 is the 
energy of quasiparticle excitations. The order parame- 
ter A and the chemical potential are fixed by the self- 
consistency equations 

A = -\J2u k v k , N = 2^v*, (7) 

k k 

where A is the BCS coupling constant. After renormal- 
ization of the coupling constant A and considering only 
the low-energy limit, the gap equation can be related to 
the two-particle scattering amplitude 0, [46| . 

The product form of the BCS wave function greatly 
simplifies the calculation of the statistics, since different 
k states can be treated separately. For a single pair of 
states (k |, —k {) the sum over all possible configurations 
can be easily performed 

e -s k (x) = (BCS|e j(fi '=T+«-'=i)x| BCS ) ( 8 ) 
= 4 + v 2 e 2 ^ . 

The sum of all states yields the result 

S( X ) = -]>>[l + ^( e 2 ^-l)] . (9) 

k 

This is one of the main results of our paper. It is valid 
in 2 and 3 dimensions; the dimension will only enter into 
the density of states, when transforming the sum over 
k into an energy integration via the standard expression 



No = m D / 2 (2e) D / 2 - 1 /2Tr D - 1 h D for D = 2,3. It should 
be noted, that in a strictly two-dimensional system the 
low-energy scattering amplitude vanishes ~ — 1 / log e and 
consequently the gap equation ([7]) shows a logarithmic 
divergence for e — » 0. However, for the more realistic 
situation of a quasi- two-dimensional cloud (i.e., a three- 
dimensional trapped atomic cloud strongly confined in 
one dimension), this singularity is eliminated. Results 
derived for the strictly two-dimensional situation are still 
valid for the quasi-two-dimensional case, however with 
the chemical potential /i shifted by the ground-state en- 
ergy 

We now discuss some limiting cases in which analytical 
expressions can be obtained. On the BEC side, /i < and 
A <c \n\ leads to v 2 <c 1 for all energies and allows to 
expand the logarithm in Eq. The result is 

S(X) = ~ (e 2lx -l) , (10) 

which corresponds to a Poissonian number statistics of 
pairs of atoms. This supports the picture of strongly 
bound pairs in a coherent state. Note, that the factor 
of 2 in the exponent leads to exponentially growing cu- 
mulants in the fermion number, viz., C n = 2 n N/2. We 
therefore expect strong fluctuations. Remarkably, the re- 
sult (|10p is in accordance with the number statistics of 
condensed bosons in a bin with volume much less that the 
total volume, see the expression given after Eq. ([SJ that 
was derived using the canonical formalism. Since we are 
counting single fermions instead of bosons in the present 
case, the counting field x is replaced by 2%, and there is a 
prefactor 1/2. This agreement is quite remarkable, since 
the starting point of our approach is the grand-canonical 
formalism. It is an indication that using the BCS wave 
function to describe the number fluctuations of a small 
subsystem works surprisingly well. 

On the BCS side the situation is quite different. Here 
fi = ep ^> A and we obtain 

S(x) = -ixN-7rND-^(\co S (x)\-l). (11) 

We observe that the first term is dominant but con- 
tributes only to the first cumulant. The fluctuations 
come from the second term in Eq. (|11[) which is smaller 
by a factor A/ep. Furthermore, similar to the degener- 
ate Fermi gas, the odd cumulants C n for n > 3 vanish, 
which is again a consequence of particle-hole symmetry. 

Due to the constant density of states in 2D we can 
obtain an analytical expression for the CGF for arbitrary 
fi and A, which reads 



S(X) = -N— 



cos(x)atan ( ^^-e lx 



. f^ep 
atan — — 



-^±ln[l + v 2 (e^-l)] 



(12) 
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Here v% = (1 + /i/vV 2 + A 2 )/ 2 is the BCS coherence 
function for k — 0. While we do not have an analyt- 
ical expression in 3D, we expect a similar behavior as 
in 2D. This will be corroborated later by comparing the 
numerically obtained cumulants in 3D to the 2D case. 




FIG. 2: (Color online) Second cumulant C2 as a function of 
£ = 1/k.Fa. For the strongly confined, quasi-two-dimensional 
system (2D), the Fermi vector can be approximated by the 
inverse of the ground state size: /cf = n/io- The dotted line 
corresponds to £ M =o- The inset shows Ci normalized to A/cf- 




FIG. 3: (Color online) Third cumulant C3 as a function of 
£ as defined in Fig. [2] The dotted line corresponds to £ M =o- 
The inset shows C3 normalized to A/ef- 

In Fig. [5] we show how the density noise changes 
through the BEC-BCS crossover. Going from the BEC 
to the BCS regime strongly suppresses the fluctuations 
in agreement with our previous discussion. The inset 
shows the fluctuations normalized to A/ep; they ap- 
proach a constant value in the BCS limit. Another im- 
portant information gained from a noise measurement 
is the order parameter, which can be extracted from 
a measurement of C2 in the BCS limit, Eq. pip , as 
A/e F = 4C 2 /ttND. Figure [3] shows the third cumu- 
lant. The global behavior is rather similar to the second 
cumulant, i.e. C3 is strongly reduced going from the 



BEC limit to the BCS limit. However, the behavior of 
C3 normalized to A/ep shown in the inset is qualita- 
tively different of that shown in the inset of Fig. [2] since 
C3 vanishes also in the BCS limit. Note also that C3 
vanishes faster in the 2D case than in the 3D case. To 
understand the behavior of the third cumulant in more 
detail we recall that C3 is related to the skewness of the 
distribution, i.e. C3 is a measure of the difference be- 
tween positive and negative fluctuations. To see this, we 
note that an elementary binomial event has the property 
ln(l + vl(exp(ix) - 1)) = i\ + M 1 + u 2 k (exp(~ix) - 1)). 
Using this property and the particle-hole symmetry of v\ 
in the BCS limit, we can rewrite the CGF in the BCS 
regime as i\N + N J deln[l + v\u\ (cos(2x) — 1)]- Here 
we have used that the density of states close to tp can 
be approximated by a constant Nq. The second term is 
even in \ and therefore only contributes to even cumu- 
lants, whereas all odd cumulants for n > 3 vanish. The 
difference between 2D and 3D seen in Fig. [3] is caused 
by the (small) energy dependence of the 3D density of 
states, which is absent in 2D. 



V. FINITE TEMPERATURES 

We would now like to discuss the effect of finite tem- 
peratures in a qualitative way. The excitations in the 
BCS theory are fermionic quasiparticles. This is a good 
approximation on the BCS side of the transition. The 
fluctuations (measured by the second cumulant C2) will 
be reduced with increasing temperature, since they ap- 
proach the linearly increasing C2 of the free Fermi gas 
that is lower at the critical temperature than C2 at 
T = 0, C f 2 lee (Tg cs ) < C 2 (0). At even higher temper- 
atures T > Tp, C2 reaches the classical Poisson value 
N. The situation on the BEC side of the crossover will 
be very different. According to our result the statis- 
tics is a Poissonian distribution of molecules and, hence, 
C2 is doubled in comparison to the classical value for 
the atomic gas [47|. Upon increasing the temperature, 
the main effect on the fluctuations will be a breaking 
up of molecules, which will reduce the second cumulant 
C2 above a dissociation temperature to finally reach the 
value for the classical gas. Thus, we expect quite a differ- 
ent temperature dependence of the number fluctuations 
on the BEC or the BCS side of the transition. 



VI. CONCLUSION 

In conclusion, we have calculated the full counting 
statistics of number densities in an ultracold fermionic 
atomic cloud with attractive interactions. The number 
statistics in the vicinity of the BEC-BCS crossover dis- 
plays interesting features which reveal the nature of the 
many-body ground state. Poissonian fluctuations of a 
molecular condensate on the BEC side are strongly sup- 
pressed on the BCS side. The size of the fluctuations in 
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the BCS limit is a direct measure of the pairing potential. 
We have also discussed the BEC-BCS crossover of the 
third cumulant and the temperature dependence of the 
second cumulant. These quantities can be accessed ex- 
perimentally and provide additional information on the 
many-body ground state in the crossover regime. The 
concept of counting statistics in ultracold gases opens in- 
teresting possibilities to study the interplay between co- 
herence and correlation in quantum many-body systems. 
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